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Abstract 

The capacity of carbon atomic chains with different terminations for hydrogen storage is studied 
using first-principles density functional theory calculations. Unlike the physisorption of H2 on 



the H-terminated chain, we show that two Li (Na) atoms each capping one end of the carbon 



chain can hold ten H2 molecules with optimal binding energies for room temperature storage. The 
hybridization of the Li 2p states with the H2 cr orbitals contributes to the H2 adsorption. However, 



the binding mechanism of the H2 molecules on Na arises only from the polarization interaction 

c . 

between the charged Na atom and the Ho. Moreover, additional H2 molecules can be bound to 
O . 

the carbon atoms at the chain ends due to the charge transfer between Li 2s2p (Na 3s) and C 

CM 

£> ■ 2p states. Importantly, dimerization of these isolated metal-capped chains does not affect the 

hydrogen binding energy significantly. In addition, a single chain can be stabilized effectively by 
G\ ■ 

termination on the Ceo clusters. With a hydrogen uptake of > 10 wt % on Li-coated C6o-C n -C6o (n 

= 5, 8), the Lii2C6o-C n -Lii2C6o complex, without reducing the number of adsorbed H2 molecules 



chains. 

03 



per Li, can serve as better building blocks of polymers than the (Lii2C6o)2 dimer. These findings 
suggest a new route to design cluster-assembled storage materials based on terminated sp carbon 
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I. INTRODUCTION 

The development of economical hydrogen energy critically depends on finding efficient 
and safe media that can store hydrogen with high gravimetric and volumetric density. 1 The 
storage materials should also reversibly adsorb/desorb H 2 under ambient thermodynamic 
conditions. These criteria not only limit the choice of storage materials to those composed 
of light elements, but also require the hydrogen adsorption energy to lie between the ph- 
ysisorbed and chemisorbed states (0.1^0.2 eV) 2 . 

Carbon nanostructures with sp 2 -like bonding-"- have long been expected to be promising 
storage materials due to their light weights and large surface areas. However, hydrogen 
molecules bind weakly to pristine carbon materials via van der Waals interactions 3 . A new 
approach, coating the carbon surfaces with early transition metals (TMs), has recently been 
shown to improve storage performance 4,5 . These complexes not only enhance the interaction 
between hydrogen molecules and TMs through the Kubas interaction (hybridization of the d 
states with the H 2 states) 9 , but also meet the U.S. Department of Energy goal of obtaining 
a gravimetric capacity of 9 wt % by the year 2015. Unfortunately, preferable clustering 
of TM atoms on the carbon nanomaterial surfaces results in poor reversibility 10 . Hence, 
the adsorption energy of metal atoms to the substrate should be larger than the cohesive 
energy of the bulk metal. The use of alkali metals (AMs) and alkaline-earth metals such as 
Li^ and Ca 12 to produce uniform coating has been proposed, as their cohesive energies are 
substantially smaller than those of TMs (~4 eV), but the binding energy of H 2 molecules 
in the Lii 2 C6o system is too small (0.075 eV/H 2 ) for room temperature application^. 

In the search for a more feasible high-capacity storage medium, we now turn to the 
5]9-hybridized carbon chains, which are considered as the important building blocks in nano- 
electronics and nano-mechanics. The reasons for choosing carbon chains in our study are 
the following: Due to their high reactivity, bare carbon chains can be stabilized by ter- 
mination on a wide variety of atoms and groups^. Isolated sp chains terminated with 
hydrogen atoms, i.e., polyynes (with alternating singlet-triplet bonds) and cumulenes (with 
double C-C bonds), have already been synthesized in their basic forms^. In addition, chains 
with metal atoms connected to the ends have been generated 15 , and their magnetic, elec- 
tronic, and transport properties have been studied extensively^^. In particular, the hybrid 
sp+sp 2 carbon-based systems 18 (with linear chains connecting sp 2 -type fragments), repre- 



senting typical interfaces in realistic nanostructures, exhibit exceptional physical and chem- 
ical properties. Recent experiments have realized single carbon chains bridging graphene 
nanoribbons 19 , and produced junctions between a single carbon chain and two single- wall 
carbon nanotubes (SWCN)^. 

In this study, we conduct theoretical studies of hydrogen storage media consisting of H, 
AMs, and Ceo terminated cumulenes or polyynes. We show that hydrogen molecules prefer 
to be physisorbed rather than chemisorbed on the H-terminated carbon chains. Thus Li (Na) 
atoms capped carbon chains are chosen to improve the hydrogen storage performance. As 
a result, H 2 molecules are adsorbed not only on the Li or Na atoms, but also on the carbon 
chains. In addition, a linear periodic polymer, composed of Li atoms and five (eighth- 
atom carbon chain bridging Ceo fullerenes, can store up to ^10 wt % of H 2 molecules 
with an average binding energy of 0.15 (0.14) eV/H 2 . Our results not only advance our 
understanding of H 2 adsorption on terminated carbon chains, but also offer a new avenue 
for efficient hydrogen storage. 

II. COMPUTATIONAL DETAILS 

Numerical calculations were performed using spin-polarized density functional theory 
(DFT) with the Perdew-Wang (1991) exchange-correlation function 21 , as implemented in 
the DMol 3 package (Accelrys Inc.). 22 A double numerical atomic orbital augmented by d- 
polarization functions (DNP) was employed as the basis set. In the self-consistent field 
calculations, the electronic-density convergence threshold was set to lxlO -6 electron/ A 3 . 
Geometric optimization was performed with convergence thresholds of 10 -5 Ha for the en- 
ergy, 2xl0 -3 Ha/ A for the force, and 10 -4 A for the atomic displacements. We performed 
normal-mode analysis of the obtained structures to ensure that the structures optimized 
without any symmetry constraints were true minima of the potential-energy surface. To 
validate the above convergence parameters, the bond length and vibrational frequency of 
the C 2 dimer were calculated to be 1.312 A and 1647 cm -1 , respectively, which agree well 
with the corresponding experimental values of 1.31 A and 1641 cm -123 . The calculated bond 
length of the H 2 molecule (0.748 A) is very close to the experimental value of 0.741 A. 



III. RESULTS AND DISCUSSION 

A. H2 adsorption on H-terminated carbon chains 

We first consider the physisorption of a hydrogen molecule on the standard cumulene 
(C5H4) and polyyne (CsH 2 ). In Fig. 1, the lowest energy configuration of a H 2 molecule is 
located above C-C bond with the H-H molecular axis between parallel and perpendicular 
to the C-C bond. The adsorption energies of H 2 to C5H4 and CsH 2 are about 0.05 eV 
and 0.08 eV, respectively, resulting from van der Walls interactions. We also calculated the 
adsorption energies of H 2 to C 2n+ iH 4 (n = 3, 4) and C 2n H 2 (n = 2, 3, 5), and found that 
the physisorption energies did not vary significantly with the chain lengths. 

Next, we investigate the H 2 dissociated chemisorption on C5H4 or CsH 2 by optimizing 
the structure with two H atoms placed close to two adjacent C atoms (i.e., CI and C2). 
Upon the atomic hydrogen adsorption, transformation of the sp hybridization of C1-C2 bond 
to the sp 2 hybridization results in a change of the bond lengths and angles. A significant 
elongation of the C1-C2 bond length in the chemisorption state indicates that this bond is 
weakened. 

We apply the nudged elastic band (NEB) method to study the minimum-energy path 
(MEP) 24 of the hydrogen molecule dissociation on the carbon chains. The image number 
considered is 16 to ensure that the obtained MEP is correct. The physisorption of the 
hydrogen molecule is chosen as the initial state (PS) and the chemisorption of two hydrogen 
atoms as the final state (CS), as depicted in Fig. 1. From PS to CS, high barriers of 3.81 eV 
and 4.32 eV must be overcome. This indicates that the dissociation process of H 2 on C5H4 
or CsH 2 is very difficult in kinetics. 

B. H 2 adsorption on Li (Na)-capped carbon chains 

The above results have demonstrated that the interaction between pristine carbon chains 
and the hydrogen molecule is very weak, such storage systems may have to be operated 
at below ambient temperature. Now the question arises: can the storage performance of 
carbon chains be improved by using alternative capping elements? We then consider several 
kinds of short atomic chains (C n =5_io) terminated with two alkali metal atoms. Our full 
structural relaxations reveal that the linear Li 2 C n and Na 2 C n structures are more stable 



and energetically favorable relative to the zigzag structures. Figure 2a illustrates that the 
average binding energies per Li (Na) atom with C n chains show an even-odd oscillatory 
behavior. These binding energies are obviously much larger than those of Li (1.63/1.52 
eV) and Na (1.11/1.41 eV) on graphene/fullerene^^. In addition, the configurations of 
Li 2 C n and Na 2 C n complexes exhibit significant differences for odd or even values of n. For 
instance, when n = 5, the C-C distances, as depicted in Figs .2b and 2d, are rather uniform. 
In contrast, the bond lengths (Figs. 2c and 2e) alternate for n = 8. This behavior implies 
that the bonding patterns of Li(Na) 2 C2 n +i and Li(Na) 2 C2 n are the same as those of C2 n +i 
and C2n, respectively. 

Note that the induced positive charge of the Li atom is smaller than that of Na, which 
is not compatible with the order of their binding energies (Fig. 2a). For example, each 
Li (Na) atom carries 0.39 (0.62) e and 0.43 (0.66) e in Li 2 C 5 (Na 2 C 5 ) and Li 2 C 8 (Na 2 C 8 ), 
respectively. We find that Li 2p orbitals become partly filled due to the carbon chain back- 
donating some electrons to Li, whereas Na still has empty 3p orbitals. This discrepancy 
means that the strong bonding between the Li atom and carbon chains is not simply ionic, 
distinctly different from the ionic bonding as the Na atom is attached to the chains. Figures 
2f and 2g clearly show that the 2p state of Li and the C p orbitals are hybridized for the 
binding of Li on the C5 and C 8 chain. Since the p-p hybridization is stronger than the s-p 
hybridization, the binding energies of Li atoms to the carbon chains are larger than those of 
Na atoms in the Na2C n complexes. This bonding mechanism has also been observed in Li 
coated boron-doped graphenes 27 . In detail, the hybridization of the Li p y orbital with the 
lowest unoccupied molecular orbitals (LUMO) of C5 gives rise to 7r-bonds (Fig. 2f), whereas 
the overlap between Li p x and C p x orbitals in Li 2 Cs forms a-bonds (Fig. 2g). Because the 
7r bond is generally weaker than the a bond, it is understandable that Li atoms bind more 
strongly to C2 n than that to C2 n +i (as revealed in Fig. 2a). 

We now turn to the interaction between the metal-capped carbon chains and hydrogen. 
The results of our calculations present that one Li (Na) atom can attach 5 H 2 molecules 
in the Li 2 C5 (Na 2 Cs) and Li 2 Cs (Na 2 Cs) complexes, suggesting that the maximum number 
of adsorbed H 2 molecules does not depend on the chain types. The optimized geometries 
for 5 H 2 molecules adsorbed on Li 2 C n and Na 2 C n (n = 5, 8), as illustrated in Figs. 3a- 
3d, demonstrate that the H-H distances are slightly enlarged compared with the free H 2 
bond length (0.74 A). All the side H 2 molecules tend to tilt toward the chain so that one 



of two H atoms of adsorbed H 2 molecules becomes relatively closer to the metal atom. To 
understand the binding mechanism of H 2 on Li- or Na-capped carbon chains, we employ 
the Mulliken charge analysis. In detail, charge transfer from the a orbital of a side H 2 to 
the attached Li atom (^0.08 e) is four times larger than that to the Na atom (^0.02 e). 
Therefore the binding of H 2 to the Na atom arises basically from the induced polarization 
of the H 2 molecules by the Na atom. However the polarization interaction is not the only 
factor responsible for H 2 binding to Li due to the average binding energies of H 2 molecules 
(Figs. 3e and 3f) in Li 2 C 5 (H 2 ) 5 (0.19 eV/H 2 ) and Li 2 C 8 (H 2 ) 5 (0.19 eV/H 2 ) being slightly 
larger than those in Na 2 C 5 (H 2 ) 5 (0.15 eV/H 2 ) and Na 2 C 8 (H 2 ) 5 (0.14 eV/H 2 ). The analysis 
of molecular-orbital coefficients for the Li 2 C5(H 2 ) 5 and Li 2 Cs(H 2 ) 5 complexes shows that 
the hybridization of the Li 2p and 2s states with the H 2 a orbitals (see Figs. 3g and 3h) 
also contribute to the H 2 adsorption. The consequences of the orbital interactions are as 
well reflected in the charge variations of the Li. The effective charge of the Li atom in both 
Li 2 C5 and Li 2 Cs complexes varies from 0.32 to 0.06 e as the number of hydrogen molecules 
increases from 1 to 5. As a result, the binding energies of H 2 molecules to the Li atom 
slightly decrease with the increasing number of H 2 molecules (see Figs. 3e and 3f). 

It is interesting to note that the H 2 molecules are adsorbed not only on the Li or Na 
atoms, but also on the carbon chains. Figure 4 shows that additional 8 H 2 molecules can be 
attached to Li 2 C5(H 2 )i , Li 2 Cs(H 2 )io, Na 2 C5(H 2 )i , and Na 2 Cs(H 2 )io complexes, resulting 
in hydrogen capacities of 32.7, 24.7, 25.4, and 20.2 wt %, respectively. The average binding 
energies of these 8 H 2 molecules to Li 2 C n (H 2 )i and Na 2 C n (H 2 )i (n = 5, 8) are 0.11 and 
0.12 eV/H 2 , respectively, which lie within the required range of 0.1-0.2 eV. The essential 
reason for hydrogen attachment to the chain is that carbon atoms participate in the electron 
transfer, and thus the carbon chain differs from the pristine one. For instance, a C atom at 
the end of Li 2 Cs (Na 2 Cs) obtains the largest amount of electrons from Li (Na) and carries 
-0.37 e (-0.46 e). Very importantly, the charge on the carbon atoms does not vary during 
hydrogen molecules adsorbed on metal atoms. Then the electric field around the C atoms, 
generated by the substantial charge redistribution, can polarize the hydrogen molecules. This 
is confirmed through the plot of the deformation electron density (Fig. 4) which exhibits 
the localized characteristics of H-H bonds and no charge transfer between H and C atoms, 
different from the charge distribution of the side H 2 molecules bound to Li atoms. The H-H 
electric dipole vector should be parallel to the electric field direction, thus the H 2 molecules 
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align vertical to the chain in such electric field. 

The above results are just promising for the isolated metal-capped carbon chains. We 
therefore put our further considerations on the cluster-assembled materials by using these 
metal-capped carbon chains as building blocks. The geometry optimization of M 2 C n =5 5 g (M 
= Li, Na) dimers was performed by starting with several initial configurations. As shown in 
Fig. 5, two M 2 C n complexes have formed stable dimers through C-C bonds or M-C bonds 
but not M-M bonds due to the repulsion between the like charges. The (Li 2 Cs) 2 dimer has a 
binding energy of 3.75 eV which is substantially smaller than 7.70 eV corresponding to the 
(Li 2 Cg) 2 dimer due to C-C bonds formation between two Li 2 Cg complexes. On the contrary, 
the binding energy of the (Na 2 Cs) 2 dimer (2.72 eV) is larger than that of the (Na 2 Cg) 2 dimer 
(1.74 eV). However the dimerization affects the number of adsorbed hydrogen molecules, as 
smaller charge on the metal atom reduces the polarization effect of the metal cation on 
the hydrogen molecules. As a result, the gravimetric densities obtained by the (Li 2 Cs) 2 , 
(Li 2 Cs) 2 , (Na 2 Cs) 2 , and (Na 2 Cg) 2 dimers are lowered to 14.0, 9.8, 13.1, and 11.3 wt %, 
respectively. Importantly, dimerizations do not change the binding energies of the H 2 on 
metal atoms significantly, and therefore should not affect the desorption temperature. 

C. H 2 adsorption on Lii 2 C6o-C n - Lii 2 C6o 

From above, we have seen that the type of chain-termination turns out to influence the 
hydrogen adsorption performance. Here we have considered the sp 2 -terminated C n chains as 
potential storage media by choosing Ceo fullerene as the end-capping candidate, because it 
not only mimics the inner SWNT caps 20 but also can hold metal atoms to attach hydrogen 
molecules. We find that C n chains, regardless of n being odd or even, energetically prefer to 
bind the common bond of two hexagons ([6,6]), resulting in the C6o-C 2n+ i-C6o and Ceo-C 2n - 
Ceo complexes with D 2 d and D 2 h symmetry, respectively. As show in Figs. 6a and 6b, the 
bond lengths of the C5 (Cg) chain in C60-C5-C60 (Ceo-Cg-Ceo) imply a polyyne (cumulene) 
structure. The binding energy is calculated by subtracting the equilibrium total energy of 
the C6o-C n -C6o complex from the sum of the total energy of isolated chain (C n ) and of Ceo 
clusters. The structure of two Ceo clusters bridged by the C5 (Cg) chain with a binding 
energy of 4.46 eV (4.94 eV) is more stable than the dumbbell-like structure of Ceo dimers 28 . 
The stability of C6o-C n -C6o (n = 5, 8) has been further verified by carrying out ab initio 



molecular dynamics simulations at 800 K with a time step of 1 fs. After running 5000 steps, 
the geometry is still kept although the chain is slightly distorted, suggesting that the Ceo- 
C n -C6o complexes are stable. This is confirmed by recent theoretical results that Ceo-Cg-Ceo 
can be stable up to a temperature of 1000 K 20 . 

Considering the Li atoms interacting with the C6o-C n -C6o (n = 5, 8) complex, it is found 
that the pentagonal face of Ceo is the preferred site. A close examination of the geometry 
of the Lii 2 C6o-C n -Lii 2 C6o configuration, as depicted in Figs. 6c and 6d, illustrates that the 
12 Li atoms adsorbed on each Ceo cluster are divided into two sets. The Mulliken charge 
analysis presents that the Li 2p orbit als are partly filled. Thus the hybridization between 
the Li 2p and C 2p orbitals strengthens the Li binding to Ceo- The average binding energy 
per Li atom with C 6 o in Lii 2 C 6 o-C5-Lii 2 C 6 o (1.74 eV) or Lii 2 C 6 o-C 8 -Lii 2 C 6 o (1.76 eV) is 
larger than that in the isolated Lii 2 C6o (1.71 eV) and the Li bulk cohesive energy (1.63 eV 
experimental value). Hence the clustering of Li atoms will not occur in our models. 

We now enter the next phase of our calculation, namely, to discuss the possible number 
of hydrogen molecules that can be bound to the Li^Ceo-Cs-Li^Ceo and Li^Ceo-Cs-Li^Ceo 
structures. Previous studies have suggested that 5 H 2 molecules can be bound to each Li 
atom in the isolated Lii 2 C6o cluster, whereas the binding energy of the fifth H 2 molecule is 
only 0.06 eV 12 . Therefore we consider if each Li atom can adsorb 4 H 2 molecules in both 
Lii 2 C6o-C5-Lii 2 C6o and Lii 2 C6o-C8-Lii 2 C6o structures, corresponding to 10.3 and 10.1 wt % 
gravimetric densities, respectively. The average binding energies of 96 H 2 molecules with 
Lii 2 C6o-C5-Lii 2 C6o (Fig- 7a) and Lii 2 C6o-C8-Lii 2 C6o (Fig. 7a) are 0.16 eV. The equilibrium 
H-H bond lengths are 0.755 A, which is comparable to that in Lii 2 C6o(H 2 ) 60 , namely, 0.753 
A. The shortest and longest distances between the Li atom and H atoms are observed to be 
of the order of 2.330-2.750 A. 

As mentioned above, a single carbon chain bridging two Lii 2 C6o clusters does not change 
the number of adsorbed hydrogen molecules per Li. In contrast, a (Lii 2 C6o)2 dimer would 
lower the hydrogen storing capacity because some of the Li atoms linking Lii 2 C6o clusters 
may not be able to attach hydrogen molecules^. Now readers have to ask whether the 
Lii 2 C6o-C n -Lii 2 C6o structures can form the building blocks of a new kind of solid storage 
materials? To address this question, we first consider the adsorption of Li atoms on linear 
periodic C6o-C n -C6o polymers. After full symmetry unrestricted geometry optimization, 
Lii 2 C6o-C n -Lii 2 C6o units can maintain their structural identity. The binding energy of Li 



atoms to a linear periodic C60-C5-C60 (Ceo-Cs-Ceo) polymer is 1.80 (1.76) eV/Li, respectively, 
indicating that the polymerization has not affected the binding strength of Li to Ceo- On 
the experimental side, one-dimensional 28 or two-dimensional 29 Ceo polymers, the hexagonal 
faces interact with each other, have been synthesized by several methods. For the fabrication 
of C6o-C n -C6o polymer, we propose a possible route to use the SWNT/single-carbon-chain 
molecular junctions^ as the initial materials. Second, each Li atom in the polymer composed 
of Lii 2 C6o-C5-Lii 2 C6o (Li^Ceo-Cs-Li^Ceo), shown in Fig. 8, can hold 4 H 2 molecules with an 
average binding energy of 0.15 eV/H 2 (0.14 eV/H 2 ), resulting in the maximum gravimetric 
density of 10 wt % (9.6 wt %). This not only meets the gravimetric density target set by 
the U.S. Department of Energy for the year 2015, but the hydrogen binding energy is also 
optimal for the system operation under ambient thermodynamic conditions. 

IV. CONCLUSIONS 

In conclusion, the adsorption of H 2 molecules on H, AMs, and Ceo terminated carbon 
chains was studied using all-electron DFT calculations. It is energetically favorable for H 2 to 
physisorb on the H terminated chains. However, for Li or Na atoms capping the chain ends, 
H 2 molecules are adsorbed not only on the metal atoms but also on the carbon atoms with 
the required adsorption energy. Although the dimerization of M 2 C n (M = Li, Na) lower the 
hydrogen uptake, the binding energy of the H 2 on metal atoms is not affected significantly. 
Thus one may assemble these metal-capped carbon chains to synthesize efficient storage 
media or catalysts. Furthermore, we propose a linear periodic polymer, consisting of Lii 2 Ceo- 
C5-Lii 2 C6o (Lii 2 C6o-C8-Lii 2 C6o), can store up to ^10 wt % of H 2 molecules with an average 
binding energy of 0.15 (0.14) eV/H 2 . We believe that current results presented here can 
be applied to other carbon chain bridging metal-coated sp 2 materials. We also hope these 
results may provide a useful reference for the design of bulk hydrogen storage materials in 
the laboratory. 
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FIG. 1: The MEPs for the dissociation process of H2 on the (a) C5H4 and (b) CsEb. The energy 
of the initial structure (PS) is set as zero. The relevant bond distances and binding energy (E#) 
are also given. 
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FIG. 2: (a) Average binding energies of two Li or Na atoms on C n as a function of the number 
of carbon atoms. The optimized bond lengths (in A) of (b) I^Cs, (c) I^Cs, (d) Na2Cs, and 
(e) Na2Cs- The large (red), medium (purple), and small (grey) balls represent the Na, Li, and C 
atoms, respectively. The panel (f) shows the Li2Cs bonding orbital results from the hybridization 
of the Li 2p orbital and C5 LUMO. (g) The bonding orbital for the Li atoms attached to Cs- The 
isovalue equals 0.015 e/A 3 . The two colors denote ± signs of the wave function. 
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FIG. 3: Optimized configurations of H2 molecules adsorbed on (a) I^Cs, (b) Na2Cs, (c) Li2Cs, 
and (d) Na2Cs along with the typical bond lengths (in A). Panels (e) and (f) show the hydrogen 
binding energies with successive additions of H2 molecules to each metal atom in M2C n (M = Li, 
Na; n — 5, 8). Panels (g) and (h) show the orbitals for the side hydrogen molecules hybridized 
with the Li 2p orbitals in Li2C5(H2)s and Li2Cs(H2)5, respectively. The isovalue equals 0.01 e/A 3 . 
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FIG. 4: Final optimized geometries of (a) Li 2 C5(H 2 )5, (b) Li 2 C8(H 2 )5, (c) Na 2 C5(H 2 ) 5 , and (d) 
Na 2 Cs(H 2 )5 holding additional eight H 2 molecules. The relevant bond distances (in A) are given. 
The right panels (a), (b), (c), and (d) show the corresponding deformation electron densities 
(molecular charge densities minus atomic charge densities). The deformed density marked in blue 
corresponds to the region that contains excess electrons, while that marked in yellow indicates 
electron loss. The isovalue equals 0.1 e/A 3 . 
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FIG. 5: Atomic configurations of dimers of (a) I^Cs, (b) I^Cs, (c) Na2Cs, and (d) Na2Cs- 
Geometries of (e) I^Cs, (f) I^Cs, (g) Na2Cs, and (h) Na2Cs dimers respectively holding 12, 12, 
12, and 18 H2 molecules. The numbers in parentheses refer to charges on metal atoms, while other 
numbers refer to bond lengths in A. 
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FIG. 6: Atomic configurations of (a) C60-C5-C60, (b) Ceo-Cs-Ceo, (c) Li^Ceo-Cs-Li^Ceoj and (d) 
Lii2C6o-C8-Lii2C6o- The numbers in parentheses denote charges on the Li atoms. The typical 
bond lengths are given in A. 
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FIG. 7: Panels (a) and (b) show the four H2 molecules adsorbed on each Li atom in Lij^Ceo-Cs- 
Lii2C6o and Li^Ceo-Cs-Li^Ceo- 
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FIG. 8: Panels (a) and (b) show the H2 molecules adsorbed on the periodic polymers composed of 
Lii2C 6 o-C5-Lii 2 C 6 o and Li^Ceo-Cs-Li^Ceo, respectively. 
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